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Abstract —Switched reluctance machines (SRMs) have 
recently become popular in the automotive market as they 
are a good alternative to the permanent magnet machines 
commonly employed for an electric powertrain. Lumped 
parameter thermal networks are usually used for thermal 
analysis of motors due to their low computational cost 
and relatively accurate results. A critical aspect to he 
modelled is the rotor-stator air-gap heat transfer, and this 
is particularly challenging in an SRM due to the salient 
pole geometry. This work presents firstly a review of the 
literature including the most relevant correlations for this 
geometry, and secondly, numerical CFD simulations of 
air-gap heat transfer for a typical configuration. A new 
correlation has heen derived: Nu = 0.181 

Index Terms — air-gap, convection, thermal model, 
switched reluctance, CFD. 

This paper is part of the ADvanced Electric Powertrain 
Technology (ADEPT) project which is an EU funded 
Marie Curie ITN project, grant number 607361. Within 
ADEPT a virtual and hardware tool are created to assist 
the design and analysis of future electric propulsion, 
especially within the context of the paradigm shift from 
fuel powered combustion engines to alternative energy 
sources (e.g. fuel cells, solar cells, and batteries) in 
vehicles like motorbikes, cars, trucks, boats, planes. The 
design of these high performance, low cost and clean 
propulsion systems has stipulated an international coop¬ 
eration of multiple disciplines such as physics, math¬ 
ematics, electrical engineering, mechanical engineering 
and specialisms like control engineering and safety. By 
cooperation of these disciplines in a structured way, the 
ADEPT program provides a virtual research lab com¬ 
munity from labs of European universities and industries 

m. 

I. Introduction 

Switched reluctance motors (SRMs) are becoming 
more popular for traction drives 121 since they are robust 
and do not require permanent magnets 1(3 . They have a 
simple rotor construction and regulated magnetic flux. 


which can provide a significantly extended speed range. 
An optimized SRM design can have power densities 
competitive with a PM machine 01. However, high 
torque ripple and associated vibration and noise are pos¬ 
sible disadvantages ||5l. In order to achieve high torque 
densities, water cooling is usually employed for traction 
motors in electric powertrains. In order to evaluate the 
cooling effectiveness, thermal analysis must be carried 
out, modelling all the paths that lead the heat from the 
active parts to the water jacket and the environment. 
Eumped parameter thermal networks (EPTN) are one 
of the most widely used approaches for thermal anal¬ 
ysis due to their low computational cost and relatively 
accurate results. One of the most critical aspects to be 
modelled is the air-gap heat transfer 0, Q, and this 
is particularly challenging in an SRM due to its salient 
pole geometry. This heat transfer path is governed by 
convection and usually the theory developed by Taylor 
ifSl is applied. Taylor studied the instability developing 
in an annulus between two smooth cylinders where the 
inner one is rotating, finding that pairs of counter-rotating 



Fig. 1. Taylor instability 















vortices appear in the z — r plane (Fig. [T]) above a 
particular rotational speed. 

As derived in Q, assuming the case of adiabatic 
flow between infinitely long smooth cylinders and very 
narrow gaps, i.e. g/Q where g is the gap width and 
Ri the inner radius of the annulus, Taylor found that the 
onset of the instability occurs at Toc = 1,697, where 
Ta is a dimensionless group called the Taylor number 
whose definition will be discussed in the next section. In 
the case of non-adiabatic flow, with heat transfer between 
the fluid flow an the boundary walls, Becker and Kaye 
10 derived and proved that the Taylor vortices should 
appear at Toc = 1, 740. After Taylor, many other authors 
explored this particular flow looking for example at the 
heat transfer between cylinders llTOl . the modes of the 
instability lITll or the effect of an additional axial flow 
m, aa, oa. a wider review of the problem may be 
found in lITSl . In general the following four modes, in 
order of increasing rotational speed, are observed Q: 

• Laminar flow 

• Laminar flow plus Taylor vortices 

• Turbulent flow plus Taylor vortices 

• Turbulent flow 

The presence of salient poles on the rotor side consid¬ 
erably changes the boundary conditions of the problem 
since the air flow is no longer driven by a continuous sur¬ 
face, and recirculations develop in the slots. This affects 
the critical value for the onset of instability and the rate 
of convective heat transfer between the rotor and stator 
surface, which may increase by up to 50% ifT^ . Thus, the 
classical correlations from the literature (e.g. II3, m, 
ifT^ l seem not to be ideal for this geometry. We focus on 
the impact of slots on the heat transfer and flow patterns 
of an SRM, by first reviewing the most relevant previous 
works and then undertaking numerical simulations of the 
airflow using the commercial CFD code STAR-CCM-i- 
v9.06 llTOl . Eventually this will enable the construction 
of a thermal resistance Rconv = representing the air- 
gap average heat transfer within a LPTN (Fig. |^, where 
A/j is the heat exchanging area and h is the heat transfer 
coefficient. Such a resistance would be a function of the 
machine geometry and speed. 

IT Literature Review 

Before starting the review we introduce the dynamic 
parameters that influence the flow field in the air gap. 
Since for vehicle traction usually an enclosed motor 
configuration is applied, no axial air flow is superim¬ 
posed to the Taylor-Couette flow in the annulus. Thus, 
the tangential Reynolds number Reg = is one of 
the main parameters that characterize the flow, where 
(jj is the rotational speed of the internal cylinder and v 



Fig. 2. Air-gap equivalent thermal resistance (averaged behaviour) 


is the kinematic viscosity of air. However, the Taylor 
number is a more appropriate dimensionless group to 
use an independent variable in annulus geometries, since 
it expresses the ratio of the centrifugal to the viscous 
forces. In our work we will use the so-called modified 
Taylor number Tam according to the definition intro¬ 
duced by Bouafia l[2Tll to take into account the presence 
of grooves: 


Uj'^Rmg^ 


( 1 ) 


Where Rm is the mean radius defined as Rm = 
h/ln{{Ri + b)/Ri) with h = Al2'KRi, and Fg is a 
geometric factor: 


TT^ Ri + Ro 
1967P 2Ri 


( 2 ) 


and 


P = 0.0571(1 - 0652^) -f 0.00056(1 - 0.652^)“^ 

rti iLi 

(3) 

The Nusselt number is the dependent variable; this 
characterizes the heat transfer by convection and the 
following definition has been used: 


Nu 


q 9 

k (Tr - Ts) 


h g 
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(4) 


where q is the specific heat flux, Tr and Ts are the 
surface temperatures of the rotor and stator side respec¬ 
tively and k is the thermal conductivity of air. We refer 
to rotor and stator sides using the subscripts R and S 
respectively. Nusselt number will be calculated as the 
average value over the entire internal surface of the 
annulus. The thermal conductivity k is taken as the mean 
of the values at the two surface temperatures. 

Throughout the literature there are a number of dif¬ 
ferent definitions for these parameters |E2l- Any system 
of definitions may be used as long as consistency is 
maintained, however for the case of slotted surfaces we 
use the system of ||2T1, as previously defined. Based on 
the various experimental geometries presented in Tab. 
we have re-expressed the correlations using a consistent 
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set of definitions of Ta and Nu in order to compare like 
with like. The values are also compared with the case 
of smooth cylinders using the correlation hy Becker and 
Kaye 111: 

Nu = 1 for Tom < Tqc 

Nu = 0.064 for Tuc < Tum < 10 ^ 

Nu = 0.205 for 10 ^ < Tum < 10 ® 

TABLE I 

Geometry Dimensions From Literature Related To Fig. [3] 

- [mm] 


Bouafla Hayase Gazley Hwang 

(2i](23i (Him 01 HI 


Case 

n.a. 

B 

c 

S1-R2 

S2-R2 

n.a. 

9 

5 

0.05 

0.05 

0.5 

0.5 

2 

Ri 

140 

0.95 

0.95 

63.3 

63.3 

34 

Lrw 

- 

- 

- 

0.1 

0.1 

c 27rRi 
^ 16 

Lrd 

- 

0.2 

- 

0.76 

0.76 

2 

OiR 

- 

10° 

- 

0° 

0° 

not given 

Nr 

0 

12 

0 

36 

36 

16 

Ro 

145 

1.0 

1.0 

63.8 

63.8 

36 

Lsw 

8.3 

- 

- 

- 

4.6 

- 

Lsd 

15 

- 

0.2 

- 

1.5 

- 

as 

0° 

- 

10° 

- 

0° 

- 

Ns 

48 

0 

12 

0 

23 

0 


It is worth noting that only Hayase 12^ . ll25l numer¬ 
ically simulated a domain geometrically similar to an 
SRM rotor. However, the simulations were performed 
considering a geometry with an higher ratio of the 
air gap to the slot height compared to a typical SRM 
configuration. Also Srinivas 1271 . l28l ran some CFD 


simulations of the SRM airgap, although assuming a 2- 
D and incompressihle domain and without reporting any 
detail about the heat transfer. 

Fig. shows that in the case of slotted surfaces in 
the laminar regime Nu is in general less than 1 and 
increasing at a small constant rate. However, keeping 
in consideration the ratios t = ^ and 6 = 
where N is the number of slots, it seems that the larger 
the slot, the lower Nu. This has also been indicated 
by Tachibana and Fukui IT^ who in their analysis also 
tested the influence of gap sizes at low rotational speeds. 

As the rotational speed increases beyond a critical 
point the heat transfer rate suddenly rises since vortices 
appear which promote fluid recirculation. However, the 
results from the literature give some uncertainty for what 
the critical Tum is. For Hwang 1261, Hayase ll25l and 
Gazley Ifl^ the instability onset seems to take place 
for Tam < Tqc- Whereas according to the data from 
Bouafla lITTIl . 1231 . the vortices appear for higher values 
of modified Taylor number, and in particular the critical 
Tum = 3,900 is measured. This discrepancy is mainly 
due to the various configurations studied which differ in 
geometry and whether the rotor is the inner or outer side 
of the annulus. 

An interesting work which has not been shown in Fig. 
1^ is Gardiner and Sabersky 12^ . where the values lie 
out of the main area of data due to the fact that the 
fluid employed was water (Prandtl number Pr = 4.5). 
However, their measurements show again a similar be¬ 
haviour of the heat transfer rate with a sudden increase 
at Tom Rs 3 X 10^. More recently the heat transfer in an 
annulus between coaxial cylinders with slots had been 
studied |[T3l . 1301 . lildl . However, these works analyse 
a Taylor-Couette-Poiseuille flow, meaning that an axial 
air flux is superimposed across the Taylor-Couette flow 
generated by the rotation of the internal or external 
cylinder. 


HI. Software validation 

Computational fluid dynamics simulations were un¬ 
dertaken using CD-adapco’s STAR-CCM-i-. First, the 
baseline scenario of heat transfer between two smooth 
cylinders with the inner one rotating and heated was 
simulated. By assuming circumferential periodicity, the 
calculations were confined to 30° in the r — 6 plane as 
assumed by Hayase (241, (251 . The axial length of the 
domain was set to allow an even number of vortices 
to develop (Ml . Additionally, to prevent end effects, a 
periodic boundary condition was applied to the axial 
ends to create an infinitely long annulus. The no-slip 
condition was applied to the outer and inner walls with 
the inner wall rotating at a constant speed. Constant 
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Fig. 4. Heat transfer vs. Taylor number for a range of previous experimental and numerical studies 


temperatures of 60 °C and 35 °C were imposed to the 
inner and outer walls respectively. In Fig. the model is 
presented for the case ofr/=^ = 0.75 and rotational 
speed large enough to let the instability grow, showing 
the velocity field characterised hy the expected counter¬ 
rotating vortices. This steady state solution was achieved 
using the k — e turbulence model and the two-layer 
all y~^ wall treatment |[20l . starting the simulation from 
the fluid at rest. As far as the fluid domain, air was 
assumed to be an ideal gas with constant specific heat 
Cair, thermal conductivity kair and dynamic viscosity 
/r = vpair, where pair is the air density. All values were 
calculated at the average temperature between the rotor 
and the stator surfaces. When it comes to simulating 
very narrow gaps, i.e. p 0.99 as in the case of the 
SRM and electrical machines in general, this turbulence 
model was unable to reproduce the instability correctly. 
As suggested in ||32l, the large eddy simulation (LES) 
approach has therefore been applied due to its ability 
to resolve smaller structures compared to the models 
based on Reynolds averaged Navier-Stokes equations. 
The WALE (wall-adapting local eddy-viscosity) subgrid 
scale model was chosen to estimate the subgrid scale 
viscosity along with the two-layer all y'^ wall treatment 
to obtain accurate results with a reduced computational 
time |[20l . The simulations were then run with this 


unsteady solver applying the same surface temperatures, 
no-slip conditions and periodic boundaries. Convergence 
was assumed to be reached when the value of A^n cal¬ 
culated as the average over the rotating surface settled to 
a constant value, which happened after 0.1 s simulation 
time. 
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Fig. 5. Taylor-Couette flow pattern rj = 0.75 

In Eig. [^we present a comparison between the numer¬ 
ical simulations using the LES model at two different 
mesh sizes, and the correlation from Becker and Kaye 
ll9l . Both the meshes were built with the STAR-CCM-i- 
using polyhedral elements: Mesh 1 consist of 262,406 
elements and Mesh 2 of 512,365; the higher number 
of elements effected only the bulk area between the 
























Fig. 6. LES model validation and mesh comparison 
TABLE II 

Srm Model Dimensions 


stator inner radius 

95.95 

mm 

rotor outer radius 

95.00 

mm 

air gap width 

0.95 

mm 

rotor pole height 

33.5 

mm 

rotor pole arc 

23° 

- 

number of slots 

6 

- 


prism layers. To correctly resolve the fields near to the 
boundaries a finer mesh consisting of 20 prism layers 
with a growing factor of 1.2 was applied, with the 
first layer height of 1.2x10“® m. This mesh resolution 
allowed always a ?/+ value helow 1. The results indicate 
that the LES model is able to reproduce the increase 
of heat transfer due to the onset of the instability. The 
results are fairly mesh independent as shown in Fig. 
which gives average Nu of the numerical solutions over 
about two rotations starting after the fluid flow settled to 
an almost steady pattern. 

IV. Results and analysis 
A. Modelling the SRM 

The case of an 8/6 SRM (8 salient poles on the 
stator side and 6 on the rotor side) was considered. The 
geometry details were taken from Q since we are mainly 
interested in SRMs for electric vehicle applications. 
These dimensions, presented in Table are in line with 
the general design rules given in ll33l . For this research 
we modelled the stator surface to be smooth (assuming 
the presence of windings and slot liner). Due to the 
periodicity of the structure, only 60° of the machine has 



Eig. 7. SRM 8/6 meshed model 


been studied, applying a periodic boundary condition 
on the tangential surfaces. Again, a periodic boundary 
condition has been applied at the axial ends to mimic 
an infinitely long annulus and prevent end effects. The 
length of the model was set to accommodate a couple of 
counter-rotating vortices. The no-slip condition was set 
on the stator and rotor surfaces which are assumed to 
be heated at a constant temperature of 60 °C and 35 °C 
respectively. To simulate the motion, the reference frame 
was set to rotate at a specific to*, the outer cylinder was 
spinning at los = —w* and the inner wall was at ujr = 0 
relative to the rotating frame. The mesh, consisting of 
1,450,394 cells, is characterized by two different scales: 
one finer, keeping fhe Mesh 1 sef up which correcfly 
simulated Taylor vortices, and one coarser for the slot 
area where larger eddies appear. As shown in Fig. [TJ 20 
prism layers was built next to the wall to ensure to 
be lower than 1. 

B. Results 

A number of rotational speeds were tested for the 
SRM domain, running each simulation using the pre¬ 
vious solution as initial condition. In Fig. the velocity 
field relative to the rotating frame is presented in the 
r — 6 at Tom = 2,802. The presence of the slot leads 
the flow to separate, generating some vortices within 
the slot. The simulations showed that in the slot a 
dominant vortex, marked with the number 1 in Fig. 



Fig. 8. Velocity field with the reference to the rotor at Tam = 2, 802 






















Fig. 9. Visualisation of Taylor vortices in the SRM airgap, Tam = 
11,208 

grows with the increasing rotor speed io. In the “2” area 
various minor vortices develop, which tend to split and 
flatten on the trailing and bottom edges of the slot with 
(jj. Although in the slot the flow has a prevailing 2D 
behaviour, in the area marked with number 3 a couple 
of counter-rotating vortices in the z — r plane were 
noticed, growing from low values of Tam (~ 300). 
The presence of recirculation flows in the slot seem to 
delay the appearance of Taylor vortices. The onset of 
the instability was not recorded at Toc but at higher 
values of Tam- Moreover, the counter rotating vortices 
were not continuous all around the circumference, but 
the simulations showed them to grow only in the gap 
between the rotor pole tip and the stator. This in turn 
affected the stability of the structures in the slot. Fig. 
shows the effect of the Taylor vortices on the temperature 
field. 

The calculated heat transfer rate is presented in Fig. 
m Due to the unsteady solver, the value of Nu is calcu¬ 
lated averaging the solution over about three revolutions 
after the transition was resolved (ss 1 rotation). We can 
see that the heat transfer for low value of Tam is below 
unity Nu and increases with u without experiencing the 
expected knee indicating a regime change. The surfaces 
bordering the slot are exposed to recirculation which 
reduces the thermal gradient and accordingly the heat 
transfer rate. As the ratio of such surfaces to the surface 
on the tip of the pole is higher than 1, the effect of the 
Taylor instability on the overall heat transfer rate is not 
noticed. 

V. Conclusions and future work 

In this work the case of heat transfer in the air gap 
of an SRM was tested. First, our CFD simulations run 
with the unsteady LES solver embedded in the com¬ 
mercial software STAR-CCM-i-, showed good agreement 
between the case of heat transfer in smooth narrow gaps. 



For the SRM configuration, the CFD results highlighted 
that with the presence of SRM rotor slots the instability 
appears for higher values of Tam, in agreement with 
what was found in previous works. It is important to 
notice that through the whole tested range of rotational 
speeds, Nu never exceeded the correlation given by 
Becker and Kaye Q. A correlation has been derived 
by least-square fitting the data: Nu = 0.181 Ta^^'^ 
for 10^ < Tam < ^ X 10^. Future work will aim at 
analysing the heat transfer at very high rotational speeds, 
i.e. 10,000 rpm which correspond to Tom oc 10^ and 
exploring the effect of different slot dimensions (e.g. 
SRM 6/4). 
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